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Abstract. We have studied the parsec-scale structure of 
the quasar 3C345 at AA1.3, 2, 3.6, and 6 cm, using to- 
tal intensity and linear polarization images from 3 epochs 
of VLBA observations made in 1995.84, 1996.41, and 
1996.81. The images show an unresolved "core" and a ra- 
dio jet extending over distances up to 25milliarcseconds. 
At 1.3 and 2 cm, images of the linearly polarized emis- 
sion reveal an alignment of the electric vector with the 
extremely curved inner section of the jet. At 6 cm, the 
jet shows strong fractional polarization (to~15%), with 
the electric vector being oriented perpendicular to the jet. 
This is consistent with the existence of a transition from 
the shock-dominated to the plasma interaction-dominated 
emission regimes in the jet, occurring at a distance of 
^1.5 milliarcseconds from the "core". 
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1. Introduction 



The QSO 3C345 (V-16; z=0.5928, Marziani et al. |1996D 
presents one of the best examples of apparent superlumi- 
nal motion in a core-dominated extragalactic radio source, 
with components traveling in the parsec-scale jet along 
curved trajectories at apparent speeds of up to 10c (Zen- 
sus et al. 1995a ). 

3C 345 has been monitored using the Very Long Base- 
line Interferometry (V LBI) t echnique since 1979 (Unwin et 



19831 Biretta et al. 1986, Brown et al. 1994, Wardle et 



|l983i Biret 
1994| , Rant; 
banov 11996 



Ranta kyro et al. 1995 , Zensus et al. 1995a , 1995b , 
Taylor \199^ , Lobanov & Zensus |l 9991 ). On 



al. 
al. 

Lobanov 

arcsecond scales the source contains a compact region at 
the base of a 4" jet that is embedded in a diffuse steep- 
spectrum halo (KoUgaard et al. 1989). 

From astrometric measurements, the parsec-scale core 
has been shown to be stationary within uncertainties of 
20/ias/yr (Bartel et al 



components) apparently ejected at different position an- 
gles (P. A. ranging from 240° to 290°) with respect to the 
jet core. The components move along curved trajectories 
that can be approximated by a simple helical geometry 
(Steffen et al. 1995 , Qian et al. 1996). The curvature of 



the trajectories may be caused by some periodic process 
at the jet origin, e.g. orbital motion in a binary black hole 
system (Lobanov |l996 ) or Kelvin-Helmholtz instabilities. 
At lower frequencies, the jet extends to the NW direction, 
turning northwards at >20 milliarcseconds (mas) distance 
from the core. 

VLBI polarimetry of nonthermal radio sources can pro- 
vide stringent information about the physical conditions 
in the parse c-scale jets. Since the early works of Cotton 
et al. ( |l984|) , substantial progress has been achieved in 
this area, particularly benefiting from the enhanced ob- 
serving capabilities of the VLBAJ^. The VLBA antennas 
have standardized feeds with low instrumental polariza- 
tion, and use high-performance receivers. A new method 
of self-calibration of th e polarimetric VLBI observations 
(Leppanen et al. 1995) has made it feasible to measure 
the fractional polarization with an accuracy of 0.15%, by 
using the target source itself as a calibrator. The tech- 
nique is non-iterative and insensitive to the structure of 
the polarization calibrator. 

In 1979-1993, 3C345 was monitored using VLBI at 
several frequencies. These observations have allowed the 
determination of the trajectories and kinematic proper- 
ties of the jet components (cf. Zensus et al. 1995a). In 



this paper we extend this effort, and present three epochs 
of multi-frequency VLBA observations which also include 
polarization information. The observations were made in 
1995.84, 1996.41, and 1996.81 at observing frequencies of 
22, 15, 8.4, and 5 GHz. We describe the observations and 
data reduction in Sect. ^. In Sect. ^, we present the total 
intensity, polarization, and spectral index images of the 
parsec-scale jet in 3C345, and discuss the properties of 
the jet emission. 

Throughout this paper, we use a Hubble constant 
-ffo=100/ikms~^ Mpc~^ and the deceleration parameter 
(70=0.5. For 3C345, this results in a linear scale of 



198^. The parsec scale jet con 3.79 /j-i pc mas~^. A proper motion of Imasyr \ then. 



sists of several prominent enhanced emission regions (jet 
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^ Very Long Baseline Array, operated by the National Radio 
Astronomy Observatory (NRAO). 
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translates into an apparent speed of /3app=19.7/i ^c. We 
the positive definition of spectral index, a {S (x v^"). 



use 



2. VLBA observations 

We carried out multi-frequency observations of 3C 345 at 
three epochs, 1995.84, 1996.41, and 1996.81, using all 10 
VLBA antennas. The observations are summarized in Ta- 
ble g. 3C 345 was observed at each frequency with a 5 min 
scan every 20 min for a total period of almost 14 h. The 
bandwidth achieved in all cases was 16 MHz. The data 
were registered in VLBA format, with an aggregate data 
rate of 64Mbits/s, recording 4 channels with 1-bit sam- 
pling (mode 64-4-1, see Romney 1985 ). Some calibrator 
scans (the calibrator sources are listed in Table % were 
inserted during the observations. 

The data were correlated at the NRAO VLBA Ar- 
ray Operations Center in Socorro, NM, US. A 2-s corre- 
lator pre-averaging time was used. Interferometer fringes 
were detected on all baselines for the target source and 
the calibrators. The data were processed using Aips|^ The 
visibility amplitudes were calibrated using antenna gains 
and system temperatures measured at each antenna. The 
fringe fitting of the residual delays and fringe rates was 
performed independently for both parallel hands of the 
correlated data, and the solutions obtained were then ref- 
erenced to a common reference antenna (Los Alamos for 
all cases). After the fringe fitting, the data were averaged 
in frequency and exported into DIFMAP (Shepherd et al. 



1994) for self-calibration and imaging. The data were fur- 



ther time-averaged, with the averaging times selected in- 
dividually at each frequency so as not to exceed the re- 
spective coherence times. We applied the clean algorithm 
and self-calibration in several cycles to produce the total 
intensity images. 

2.1. Data analysis of the total intensity data 

In order to quantify the source structure, we fitted models 
consisting of elliptical Gaussian components to each of 
the self-calibrated data sets, using the tasks modelfit 
and MODFIT of the Caltech Package (Pearson & Readhead 



198S ) . For the purpose of studying the spectral properties 



of the jet, we reproduced all images in a consistent fashion, 
by selecting matching itw-ranges at all frequencies and by 
restoring the final images with the same circular beam. 

2.2. Imaging the polarized emission 

Images of the linearly polarized intensity^ were obtained. 
The instrumental polarization was determined using the 
feed solution algorithm developed by Leppanen et al. 



^ Astronomical Image Processing System, developed and 
maintained by the NRAO. 



P = Q + iU = pe^''^ = mle'''^, where Q and U are the 
Stokes linear polarization parameters, p = ml is the polarized 
intensity, m is the fractional linear polarization, and x is the 
position angle of the electric vector (EVPA) in the sky. 



( 1995 ) . The average accuracy of the antenna D-terms was 
better than 1%. For some antennas, larger values of the D- 
terms were found at higher frequencies: Brewster — 4% at 
22 GHz; OVRO— 6-8% at 22 GHz and up to 4% at 15 GHz; 
Saint Croix— 5% in 1996.81 at 22 GHz. 

We transferred the clean component models from the 
total intensity (/) images back into AiPS and then used 
the task IMAGR to produce the Stokes U and Q images 
of 3C345. The distribution of the polarized flux p and 
the position angle of the polarized emission was obtained 
from the combination of the U and Q images. The reg- 
istration of the / and P maps to within a small fraction 
of a fringe spacing was assured by using the same set of 
antenna phases in the U and Q imaging. 

The final step in the polarization calibration was to de- 
termine the absolute offset between the right and left hand 
phases at the reference station. This can be done in one of 
three possible ways: a) by observing a radio source with 
a known electric vector position angle x, b) by comparing 
the VLBI maps with VLA maps (that should account for 
the total polarized intensity), and/or c) by comparing the 
values of leakage of RCP into the left feed and LCP into 
the right feed in the receivers (quantified as "D-terms" ) of 
the antennas with values from neighbor epochs (relying in 
their time stability). We have applied the approaches b) 
& c) in our data reduction. For the epochs and frequencies 
without such information, the right- and left- hand offsets 
were introduced ad hoc, based on the consistency with the 
calibration at other frequencies and epochs (for instance, 
we were assuming that the outer, extended components 
cannot change their x very rapidly from one epoch to an- 
other). We summarize the polarization calibration for all 
epochs and frequencies in Table |[ 



Table 2. Polarization calibration summary 

22 GHz 15 GHz 8.4 GHz 5 GHz 

1995.84 D-term'^/VLA'" D-term.^/VLA'^ VLA'^ VLA" 
1996.41 VLA'' Man.'= VLA'^ Man.= 

1996.81 D-term'^ D-term'' Man.'^ Man.'^ 

^ Calibrated comparing with the D-terms at epoch 1995.83 
(Alberdi, priv. comm.). 
Calibrated using VLA images. 

The right/left hand offsets were introduced manually, based 
on the consistency with the calibration at other frequencies 
and epochs. 

Calibrated comparing with the D-Terms at epoch 1996.85 
(Gomez, priv. comm.). 



At frequencies >5 GHz, the ionospheric effects are neg- 
ligible and can be safely ignored. We also do not ap- 
ply Faraday rotation corrections to the observed electric 
vector position angles, relying on the results of Taylor 
( 1998 ), who reports small rotation measurements {RM = 
— 130radm~^ at the core and -70radm~^ at the parsec- 
scale jet, derived from a least-squares fit to the polariza- 
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Table 1. 3C345 observation and map parameters 

3C 345 image parameters 





Freq. 






Beam 












Epoch 


V 


Calibrators 


Beam size 


P.A. 


Q a 
Otot 


S \^ 


Contours in 


Fig. 1 
t> □ 




[GHz] 




[mas] 




[■ly] 


[Jy/boam] 


[(mJy/boam)x 


(lovols)] 


1995.84 


22 


3C279 


0.480x0.326 


-15?0 


6.98±0.14 


2.56±0.05 


5.1x 


(-1,1,1.4,- ■ 


.,256,362) 




15 


3C279 


0.690x0.453 


-0?7 


7.98±0.16 


3.66±0.07 


7.2x 


(-1,1,1.4,- ■ 


.,256,362) 




8.4 


3C84 


1.254x0.731 


2?9 


9.32±0.18 


4.46±0.09 


6.7x 


(-l,l,1.4,■■ 


.,362,512) 




5 


3C 84 


2.443x1.348 


-21?0 


7.51±0.15 


4.41±0.09 


4.4x 


(-1,1,1.4,.. 


.,512,724.1) 


1996.41 


22 


NRAO 91 


0.433x0.316 


-3?1 


4.51±0.09 


1.80±0.04 


3.7x 


(-1,1,1.4,.. 


.,256,362) 




15 


NRA091 


0.624x0.457 


-7? 3 


6.70±0.13 


2.67±0.05 


5.3x 


(-1,1,1.4,.. 


.,256,362) 




8.4 


OQ208 


1.136x0.810 


-7? 9 


7.49±0.15 


3.55±0.07 


5.3x 


(-1,1,1.4,.. 


.,362,512) 




5 


OQ208 


1.967x1.318 


-9? 5 


7.21±0.14 


3.91±0.08 


3.9x 


(-1,1,1.4,.. 


.,512,724.1) 


1996.81 


22 




0.451x0.322 


-15?5 


4.61±0.09 


2.09±0.04 


4.2x 


(-1,1,1.4,.. 


.,256,362) 




15 


3C286 


0.614x0.474 


-6?1 


5.98±0.12 


2.72±0.05 


5.4x 


(-1,1,1.4,.. 


.,256,362) 




8.4 


3C286 


1.132x0.887 


-16?8 


7.52±0.15 


3.59±0.07 


5.4x 


(-1,1,1.4,.. 


.,362,512) 




5 


3C286 


1.861x1.532 


-15?3 


7.26±0.15 


4.29±0.09 


4.3x 


(-1,1,1.4,.. 


.,512,724.1) 



" Total flux density recovered in the map model. The calibration errors are always smaller than a 2%, which translates into the 

tabulated values for the uncertainties in 5tot and 5peak. 
^ Peak of brightness in the map. 



tion angle measurements between 8.1 and 15.2 GHz) at 
frequencies higher than 5 GHz. 

3. Results and discussion 

3.1. Total intensity images and model fits 

The total intensity images of 3C 345 are shown in Fig. |l|. 
The parameters of the images are given in Table The 
quality of the data and the performance of the VLBA 
antennas and correlator were excellent. The model fits for 
all images are presented in Table ^. 

The errors in component position determined at each 
epoch are typically smaller than a tenth of a beamwidth, 
but to be conservative we have assumed a magnitude 
of a fifth of beamwidth. That means, approximatively, 
0.08 mas for the 22 GHz results, 0.11 mas for 15 GHz, 
0.19 mas for 8.4 GHz, and 0.32 mas for 5 GHz. The er- 
rors for the other parameters were determined combining 
the statistical standard errors provided by the task er- 
FIT (revised and improved version of errfit) from the 
Caltech Package, and some statistical and empirical con- 
siderations. 

3.1.1. The models at 22 GHz 

Table ^ lists the Gaussian fits for the data at 22 GHz. 
We can model the radio source reliably with four compo- 
nents, that we associate with the core D, the components 
C8 and C7, and a more extended emission to the NW la- 
beled as "jet" (which contains the components identified 
with C5 and C4 at lower frequencies). Components C8 
and C7 move away from the core, at angular speeds of 
0.27±0.12mas/yr (estimated from our three observations, 
and corresponding to an apparent speed of 5.3±2.3/i~^ c, 
uncertainties shown here and in the next subsections are 
statistical standard errors from a weighted linear fit to 



the core distances) and 0.35±0.12mas/yr (7.0±2.3/i^^ c), 
respectively. In Fig. ^ we represent their proper motions 
by the changes of r (distance to component D) with time. 
The inner components of the jet of 3C 345 have displayed 
complex trajectories in previous epochs. We cannot asso- 
ciate C8 with previous components, since it is new, com- 
ponent C7 has been observed before and displays acceler- 
ation from 1991 to our observing epochs. With our time 
span of a year we linearize the behavior of the C7 and 
C8 components to compute their apparent velocities. The 
comparison of the positions reported here with later ob- 
serving epochs at the same frequencies should constrain 
the parameters of their kinematics. 
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Fig. 2. Core separations for the jet components C7, C8, 
and C9 obtained from the model fits to the data at 15 and 
22 GHz. 
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Table 3. Model fitting results 







5^ 




6" 














[mas] 


n 


[mas] 


\°] 


22 GHz 


1995.84 


Core (D) 


2.730±0.140 






0.13±0.01 


1.0 






C8 


2.252±0.110 


0.40±0.08 


-122±11 


0.25±0.01 


0.66±0.01 


47±1 




C7 


1.214±0.060 


1.20±0.08 


-95±4 


0.51±0.01 


0.55±0.01 


71±1 




Jet 


0.400±0.020 


6.20±0.14 


-68±1 


4.03±0.14 


0.61±0.03 


58±3 


1996.41 


Core (D) 


2.168±0.110 






0.16±0.01 


1.0 






C8 


1.184±0.060 


0.48±0.08 


-118±9 


0.30±0.01 


0.60±0.01 


-23±1 




C7 


0.884±0.050 


1.36±0.08 


-93±3 


0.53±0.01 


0.60±0.01 


66±1 




Jet 


0.334±0.017 


5.21±0.11 


-73±2 


6.8±0.2 


0.33±0.02 


-64±1 


1996.81 


Core (D) 


1.618±0.200 






0.14±0.01 


1.00 






C9 


1.082±0.200 


0.16±0.08 


-100±26 


0.07±0.01 


0.5±0.1 


119±8 




C8 


1.087±0.050 


0.67±0.08 


-108±7 


0.39±0.01 


0.67±0.01 


122±1 




C7 


0.694±0.040 


1.55±0.08 


-93±3 


0.70±0.01 


0.52±0.01 


59±1 




Jet 


0.230±0.012 


5.9±0.1 


-72±4 


5.0±0.2 


0.5±0.2 


-51±2 


15 GHz 


1995.84 


Core (D) 


2.616±0.200 






0.15±0.01 


1.0 






C8 


3.135±0.200 


0.38±0.11 


-120±16 


0.28±0.01 


0.45±0.08 


-129±1 




C7 


1.779±0.090 


1.19±0.11 


-95±5 


0.59±0.02 


0.46±0.01 


74±1 




Jet 


0.556±0.028 


5.8±0.3 


-74±2 


6.2±0.1 


0.30±0.01 


-56±1 


1996.41 


Core (D) 


2.501±0.130 






0.12±0.01 


1.0 






C8 


2.227±0.110 


0.46±0.11 


-115±13 


0.32±0.01 


0.85±0.01 


7±1 




C7 


1.544±0.080 


1.36±0.11 


-93±5 


0.60±0.01 


0.58±0.01 


66±1 




Jet 


0.440±0.022 


6.0±0.3 


-72±2 


4.0±1.0 


0.5±0.2 


128±1 


1996.81 


Core (D+C9) 


2.844±0.140 






0.16±0.01 


1.0 






C8 


1.762±0.090 


0.57±0.11 


-108±11 


0.38±0.01 


0.79±0.01 


97±1 




C7 


1.017±0.050 


1.48±0.11 


-93±4 


0.70±0.01 


0.50±0.01 


59±1 




Jet 


0.374±0.020 


5.89±0.6 


-73±3 


4.9±0.1 


0.39±0.01 


-55±1 


8.4 GHz 


1995.84 


Core (D) 


2.746±0.197 






0.22±0.10 


1.0 






C8 


2.337±0.185 


0.33±0.19 


-119±30 


0.21±0.02 


0.72±0.08 


106±22 




C7 


2.490±0.195 


1.09±0.19 


-93±9 


0.59±0.01 


0.57±0.01 


66±1 




C5b 


0.082±0.020 


2.35±0.3 


-92±10 


1.0±0.2 


0.5±0.2 


100±9 




C5a 


0.141±0.010 


3.76±0.19 


-79±3 


1.16±0.07 


0.77±0.07 


105±11 




C4 


0.552±0.028 


6.29±0.3 


-69±7 


3.0±0.2 


0.7±0.1 


132±2 




CI 


0.203±0.010 


11.0±0.3 


-72±9 


8.4±1.0 


0.5±0.4 


112±4 


1996.41 


Core (D) 


2.097±0.205 






0.24±0.01 


1.0 






C8 


2.304±0.115 


0.43±0.19 


-115±24 


0.32±0.01 


0.59±0.06 


-27±2 




C7 


2.081±0.105 


1.28±0.19 


-93±8 


0.63±0.01 


0.57±0.01 


67±1 




C5 


0.124±0.010 


2.26±0.19 


-93±5 


1.5±0.5 


0.5±0.3 


117±2 




C4 


0.713±0.036 


6.07±1.0 


-71±8 


4.36±0.03 


0.50±0.01 


126±1 




CI 


0.108±0.010 


11.8±0.19 


-73±6 


6.3±0.2 


0.41±0.03 


141±2 


1996.81 


Core (D) 


2.154±0.110 






0.15±0.01 


1.0 






C8 


2.870±0.130 


0.54±0.19 


-107±20 


0.4±0.2 


0.84±0.01 


127±20 




C7 


0.893±0.050 


1.31±0.19 


-90±9 


0.51±0.01 


0.60±0.02 


63±2 




C5b 


0.767±0.038 


1.66±0.19 


-97±7 


0.7±0.2 


0.76±0.01 


49±4 




C5a 


0.171±0.015 


3.75±0.19 


-84±8 


2.45±0.04 


0.5±0.1 


106±1 




C4 


0.556±0.030 


6.60±0.19 


-69±6 


3.4±0.5 


0.66±0.01 


119±1 




CI 


0.131±0.010 


12.17±1.0 


-72±7 


7.0±2.0 


0.6±0.2 


141±5 


5 GHz 


1995.84 


Core (D) 


1.892±0.345 






0.15±0.10 


1.0 






C8 


1.695±0.328 


0.37±0.32 


-112±41 


0.49±0.01 


0.4±0.2 


87±3 




C7 


2.394±0.130 


1.08±0.32 


-90±17 


0.63±0.01 


0.4±0.2 


64±1 




C5 


0.189±0.010 


3.20±0.32 


-82±6 


2.2±0.1 


0.4±0.2 


97±4 




C4 


0.800±0.040 


6.30±0.32 


-69±3 


3.58±0.03 


0.54±0.01 


120±1 




CI 


0.305±0.020 


17.6±0.5 


-59±2 


17.2±0.6 


0.26±0.02 


152±2 


1996.41 


Core (D) 


1.100±0.153 






0.15±0.10 


1.0 






C8 


2.279±0.150 


0.39±0.32 


-115±40 


0.3±0.1 


0.5±0.1 


-68±19 




C7 


2.343±0.150 


1.26±0.32 


-93±14 


0.69±0.05 


0.43±0.03 


68±3 




C5 


0.420±0.030 


3.76±0.32 


-85±5 


3.8±0.2 


0.37±0.01 


108±1 




C4 


0.700±0.040 


6.65±0.32 


-68±3 


3.50±0.03 


0.60±0.01 


104±1 




CI 


0.356±0.010 


17.3±0.4 


-61±2 


15.0±1.0 


0.49±0.02 


142±2 


1996.81 


Core (D) 


0.948±0.100 






0.15±0.03 


1.0 






C8 


2.929±0.150 


0.49±0.32 


-112±33 


0.49±0.01 


0.66±0.01 


93±1 




C7 


2.025±0.100 


1.42±0.32 


-94±13 


0.84±0.01 


0.47±0.01 


64±1 




C5 


0.347±0.020 


4.03±0.32 


-82±5 


3.44±0.02 


0.4±0.2 


104±1 




C4 


0.677±0.040 


6.79±0.32 


-68±3 


3.57±0.02 


0.49±0.01 


104±1 




CI 


0.368±0.020 


16.7±0.5 


-62±2 


17.7±0.3 


0.36±0.01 


141±1 



" S, flux density of the Gaussian component. ^ r, distance to the main component. 6, position angle (doflnod north to west) 
of the component. a, major a^xis at the FWHM of the elhptical Gaussian. " b, minor axis; b/a ratio. ' (j), position angle of the 
major axis. 
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Fig. 1. VLBA total intensity flux images of 3C 345 for the three epochs and the four frequencies. The map parameters 
are given on Table 1. The synthesized interferometric beams are represented at the bottom left of each image. Notice 
that the scales are not the same at different frequencies. 



A better fit is obtained at the epoch 1996.81, by in- 
troducing a new component between the core D and C8, 
labeled as C9, very close to D, at 0.16 mas (P. A. ~ -100°). 
The model fitted fluxes at D and C9 are strongly corre- 
lated. Different values for the flux densities of D and C9 
(preserving constant the sum of both) provide almost the 
same values for the agreement factor in the fit. To the first 
order of accuracy, the fluxes at the core and C9 are close 
to 5d ~ 1-6 Jy and Scg ~ 1.0 Jy. 



We show the flux density evolution of the different 
components in Fig. ^ comparing also the values with the 
total flux density from the maps (Table |l|) and the sin- 
gle di sh me asurements obtained in Metsahovi (Tera sranta 
et al. 199J , and priv. comm. ). Lobanov & Zensus ( |1999| ) 
have applied a flare model to the observed variations of the 
22 GHz flux density of the core during our observations. In 
their description, the core is at the late stages of the decay 
after a flare in 1995.2. The new component C9 may there- 
fore be related to the flare in 1995.2. Similarly, previous 
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components have been associated to flares. The periods 
of emergence of components are r elated to the periods of 
variability (see, e.g., Zhang et al. [l998[) . T his b ehavior is 
extremely clear for 3C273 (Tiirle r et a l. |l999|) . Follow- 
ing the work of Villata & Raiteri ( |1999D on MknSOl, the 
most enhanced emission regions in the jets (called "com- 
ponents") would be the result of the boosting of helical 
features caused by the viewing angle. The helicity would 
be generated by the rotation of a binary black hole sys- 
tem. A similar model was discussed for QSO 1928-f 738 by 
Hummel et al. (1992). The periodicity in the helix would 
then be correlated with the light-curve periodicity in the 
optical and radio, in this scenario. 
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Fig. 3. Evolution of the flux densities at 22 (Terasranta 
et al. 1998D and 15 GHz (Univ. of Michigan) for the single 
dish flux, the total VLBl flux (from Table [T]) and the 
fluxes of the model flt components (from Table | 



The labeling of components follows Lobanov ( 19961 ). 
To identify our components with the trajectories reported 
before, we have combined our component positions with 
previous ones. We show in Fig. ^ our data combined with 
those mentioned above. There we can track C7, see how 
C6 is not detected now, and also the presence of C8 and 
C9. The inspection of the flux density values for C7 in pre- 
vious data and our data shows a continuity in the decay. 



from 6 Jy around 1992 to the values presented above. The 
identiflcation, thus, is unambiguous. 
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Fig. 4. Core separations of the jet components C6, C7, C8, 
and C9 at 22 GHz obtained from our model flts (last three 
epoch s) and previous results (Lobanov 1996 and Leppanen 
et al. 1995). The evolution and acceleration of C7 is ob- 
vious. The component C6 is not assigned in our model 
fitting, and C8 and C9 are new with respect to the previ- 
ous data sets. 



3.1.2. The models at 15 GHz 

The model fitting of the source structure at 15 GHz follows 
the same procedure, and provides results similar to the 
models at 22 GHz. Here, we also identify the core compo- 
nent D, two jet components C8 and C7, and the extended 
component related to the faint emission of the outer com- 
ponents of the jet. We show the results in Table |[ 

The introduction of a new component C9 at epoch 
1996.81 between the core D and C8 does not improve the 
fit. We plot the flux density changes in the lower panel of 
Fig. ^ comparing them with the total mapped flux den- 
sity and single dish measurements from the University of 
Michigan Radio Astronomy Obser vatory (hereafter UM- 
RAO, see, e.g., Aller & AUer |l996| ). 

The proper motions of C8 and C7 measured at 15 GHz 
(plotted in right panel of Fig. ||) are 0.19±0.16mas/yr 
(3.8±3.1/i-ic) and 0.30±0.16mas/yr (5.9±3.1 /I'^c), re- 
spectively. 



3.1.3. The models at 8.4 GHz 

Model fitting of the source structure observed at 8.4 GHz 
(and also at 5 GHz) requires the inclusion of the inner com- 
ponents seen at the higher frequencies, and additionally 
of several more extended components at larger distances 
from the core. We present these resulting model fits in 
Table |[ Apart from the core D, C8 and C7, a strong com- 
ponent (~0.6 Jy) is found at ~6.2mas (P. A. ^ —70°) away 
of the core for the three epochs. It is most plausibly iden- 
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t ified with the C4 component discussed by, e.g., Lobanov 
(119961) . The models are not robust for the region between 
C7 and C4, and provide one or two components for dif- 
ferent epochs. For consistency with other frequencies, we 
labeled it as C5 or C5a and C5b (if double). 

A new component is needed to fit also the emis- 
sion at larger distances. The extended emission to 
the NW at >10mas can be also reproduced with 
an extended component (FWHM of ~7mas), labeled 
as CI. The extended nature of these jet components 
makes proper motion determinations very uncertain. 
Tentative values would be, for C8, 0.21±0.28mas/yr 
(4.2±5.4/i-ic), C7: 0.23±0.28mas/yr (4.6±5.4 /i-^c), and 
C4 with 0.33±0.37mas/yr (6.5±7.2 /i^^c). We display 
these proper motions in Fig. |^. 
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Fig. 5. Distance to core of components C8, C7, C5, and C4 
at 8.4 and 5 GHz from our model fit. C5 splits into C5a 
and C5b in two epochs at 8.4 GHz, where the modeling 
does not provide robust fits to discern between one or two 
components. 



3.1.4. The models at 5 GHz 

We modeled the visibility data at 5 GHz with the same 
set of Gaussian components as in the models for 8.4 GHz. 
The results are presented in Table ^. The components 
in the core region D/C8/C7 are strongly correlated in 
flux density with each other, which limits the accuracy 
of the flux densities. However, their positions can be de- 
rived accurately. The component C4 is evident at 
(P.A. ^ —70°). The emission between the core region 
and C4 can be reproduced with a component at about 
3.5 mas out of the core. Also, the most extended emis- 
sion of the jet, is reproduced with the component la- 
beled CI, located at -15 mas (P.A. 60°). Wc show 

the evolution of the core separations of C8, C7, C5 and 
C4 in Fig. |. Values for the proper motions are C8: 
0.12±0.46mas/yr (2.3±9. 1 /i-^c), C7: 0.35±0.46mas/yr 
(6.8±9.1 /i-ic), C5: 0.86±0.46mas/yr {17. 0±9. 1 h-^c), 
and C4: 0.51±0.46mas/yr (10.0±9.1 h-'^c). (Note that the 



uncertainties presented are formal statistical standard de- 
viations from a weighted linear fit to the three points for 
every component.) 

3.2. Radio spectra 



Lobanov (1998b) showed that it is possible to map the 
synchrotron turnover frequency distribution using nearly 
simultaneous, multi- frequency VLBI observations. The 
turnover frequency is defined as the location of maxima 
{Sm,i^m), in the spectrum of the synchrotron radiation 
from a relativistic jet (spectral shape Sl^) cx v"). The 
spectrum is characterized by this maximum and the two 
spectral i ndices, athick (j^ ^ I'm) and athin {i^ > t'm)- 
Lobanov (19981) showed the first results from the imag- 
ing of the turnover frequency distribution of 3C 345 from 
observations at 1995.48, at the same frequencies as those 
reported in the present paper. 

We produced special renditions of our maps at all fre- 
quencies with a circular convolving beam of 1.2 mas in 
diameter (see Fig. ^), in order to produce spectral index 
maps of 3C345 (Fig. 0). We display only the images for 
epoch 1996.41. In Fig. ^ we also show cuts through the 
spectral index maps along a line with a P.A. of 117°. 
Epochs 1995.84 and 1996.81, which are not represented 
for brevity, show very similar features. 

The registration of the images was done by consider- 
ing the relative separations of the model fit components 
from the core (assuming that the components are optically 
thin), the theoretical predictions for the core registration 
(see Lobanov |1998a ), and the alignment of the jet compo- 
nents directly in the images. 

In the spectral index maps, two different regions in 
the parsec-scale structure of 3C345 are evident. The 
D/C8/C7 region is optically thick. The region of the more 
extended jet from about > 1.5 mas out of the core to the 
NW (P.A. ~ -70°) is optically thin. The transition be- 
tween both areas is located to the west of the component 
C7. 

The core region displays spectral indices that allow one 
to deduce that the turnover frequency at this region is 
above 15 GHz, peaking between the core and C8, and de- 
creasing to the west. The component trajectories obtained 
exhibit apparent proper motions of 0.2-0.3 mas/yr. The 
peaks of the spectral images are consistent with a shock 
evolution. This behavior and the shock parameters have 
been studied in detail by Lobanov & Zensus (1999) for 
previous observing epochs. 

In the comparison of data sets of 5/8.4 and 15/22 GHz 
we find a gradient of the spectral index, increasing to the 
south: the southern edge of the jet is optically thicker than 
the northern edge. This has been observed at all three 
epochs, and we reject the possibility of a registration error 
after a careful test of different shift values. In the maps at 
22 GHz, we observe a steeper fall in the flux density from 
the components to the south than to the north. There is 
some extended emission in this northern region. These gra- 
dients can be interpreted as an effect arising from changes 
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Fig. 6. Special renditions of the 
maps of 3C 345 for epoch 1996.41 
with a circular restoring beam 
of 1.2 mas diameter used to pro- 
duce the spectral index maps 
shown in figure |^. The contours 
are 3 mJyx (-1,1,1.73,3,...,). The 
peaks of brightness are of 3.69, 
4.24, 4.30, and 3.05 Jy/beam re- 
spectively at 5, 8.4, 15, and 
22 GHz. The shift of (10,-10 mas) 
with respect to the maps of 
Fig. |l] was made to reduce map 
sizes by computing reasons (the 
shift is unimportant as relative 
right ascension and declination 
are shown). Epochs 1995.84 and 
1996.81 show similar features. 
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Fig. 7. Spectral index maps (top in every subimage) and profiles (bottom) among the maps convolved with a 1.2 mas 
circular restoring beam (Fig. g) from our multifrequency data set at 1996.41. The profile line follows the diagonal line 
over the maps, with a P.A. of 117°. 



in the speeds across the jet flow. The curved trajectories 
of the components can also be related to this effect. 

At larger distances along the parsec-scale jet 
(>1.5mas), the spectrum shows an optically thin region, 
with the spectral index decreasing, becoming as steep as 
^22 GHz ^2. The peaks in the maps that include 5 GHz 
information can also be affected by beam overresolution. 
A similar gradient to the south is also evident in this re- 
gion in most of the cases. There seems to be an intrinsic 
asymmetry in the jet, as we will also see in the polar- 



ization analysis, in Sect. 3.3. For a helical jet geometry 
(e.g., Villata & Raiteri 199£) this intrinsic asymmetry can 
be explained by the beaming of emission towards the ob- 
server for one side of the helix — the southern in our case, 
and fainter emission in the other side. 

In summary, the spectral behavior of the jet shows 
a scenario with a clear transition between two zones with 



different physical properties: the core region and the outer 
parsec-scale jet region. The components in the latter show 
a slower flux decrease than in the inner region (comparing 
the values in the model fit), and also higher speeds of the 
component proper motions. 



3.3. Polarimetric imaging 

Before discussing the VLBI results, we want to point out 
some facts observed in the single-dish monitoring results 
of the polarized emission of 3C345. In Fig. ||, we plot 
the degree of polarization m and the EVPA x from the 
UMRAO data at epochs close to our observations. 

The polarization angle of the total emission x shows 
a value of ~ 40° close to 1995.84 and is constant around 
~ 20° after early 1996. At higher frequencies, x changes, 
showing a decrease and a further increase in its value. 
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Fig. 8. Single dish m and x of 3C345 (UMRAO). 

Xi5GH^° changes from values of - 10° at 1995.8 to 40° 

in 1996.6, turning back to ~ 20° on early 1997. The degree 
of polarization remains between 1-2% during our three 
observing epochs. Interestingly, it increases up to 4% after 
our last VLBI observations in 1997.1, when the values x 
at the three frequencies coincide again. This is the epoch 
where the C9 component appears in our model fitting, and 
also when the flux density begins to increase at the higher 
frequencies. 

We can assume that the single dish polarization mea- 
surements of 3C 345 are almost equivalent to integrating 
all of the parsec-scale jet polarized components. The VLA 
polarized images mentioned above are consistent also with 
the single-dish measurements. At the epoch of 1995.8, 
X^?GHz-5°,and xy^^c^Hz -10°. 

In Fig. ^, we show the composition of the total inten- 
sity / (contours), the polarized intensity p (grey scale) 
and the electric vector orientation angle x (segments, 
length proportional to p) for the observations in 1995.84 
at 22 GHz. It is obvious that the electric vector is aligned 
with the extremely curved jet direction in the inner 3 mas. 

Fig. ^ shows the total and polarized flux density for 
the three epochs along two cuts through the core, which 
are oriented at P. A. of 60° and 90°, which cross, respec- 
tively, the components C8 and C7. The evolution of the 
total flux for the core D and components C8 and C7 is 
evident in the figure. The brightness peak of the polarized 
emission moves westwards, and a new component appears 
in 1996.81. The evolution of C8 is the most dramatic: in 
1995.84, it dominates the p-map, then it becomes much 
weaker in 1996.41 and vanishes in 1996.81. The polariza- 
tion angle of C8 rotates from the 1st to the 2nd epoch. 
The component C7 shows a westward motion, similar in 
total and polarized emission. The value of x for this com- 
ponent is nearly constant throughout the period of obser- 
vations. We can compare those results with the 22 GHz 
map fr om 3C 345 at epoch 1994.45 presented in Leppanen 
et al. (1995), 1.39yr before our first epoch. Component 



C8 had emerged and was clearly visible in the polariza- 
tion map, but not in the total intensity image. Its EVPA 
began to be oblique with respect to the values of D and 
C7, which were parallel to the jet, as it is observed also 
in our data. Our maps at 22 GHz represent a later stage 
in the evolution of the C8 component, which is rotating 
through its evolution. 
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Fig. 10. Logarithmic slice profile of the total and polar- 
ized flux density maps of 3C 345 at 22 GHz through the 
core and the C8 (top) and C7 (bottom) components, cross- 
ing the center of the maps with a P.A. 60°and 90°, respec- 
tively (shown at the top left of each flgure). The motion of 
the components in total intensity can be appreciated, and 
also the changes in the polarized flux density features. 



Gomez et al. ( |1994[ ) have modeled features similar to 
C8 from simulations of the emission in a bent shocked rel- 
ativistic jet, reporting anticorrelation between the polar- 
ized and the total intensity emission, which is observed at 
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Fig. 9. VLBA /, p and x 
images of 3C345 at 22 GHz, 
epoch 1995.84. The total inten- 
sity / is represented with con- 
tours (value of 6mJy/beamx — 
1,1,2.24,5,11.18,25,...), super- 
imposed over a grey scale po- 
larized intensity map (peak of 
brightness of 112.4mJy) and the 
superimposed electric vectors (x, 
length proportional to p, 1 mas 
in the map is equivalent to 
lOOmJy/bcam). 



early stages of evolution of a relativistic shock. We should 
point out that a possible superposition between compo- 
nents with different electric vector orientations may lead 
to cancellation of the polarized flux and can produce ap- 
parent gaps between polarized emission components. The 
brighter polarized emission in C8 might be explained by a 
shock wave in a curved jet. More generally, the features in 
Fig. ^ ca n be e xplained in terms of the shock model (War- 
die et al. 1994) in the framework of a helical geom etry for 
the motio n of t he components (Stef fen et al. 1995| , Villata 
& Raiteri 1999). Lister et al. (1998) report similar vector 
alignment with the jet at 43 GHz, for the inner jet regions 
in several blazars. 

In Fig. |ll|, we display the images of polarized emission 
in 3C345 at 15 and 8.4 GHz, following the same scheme 
as in Fig. ^. The same features seen at 22 GHz are present 
here. At 8.4 GHz, weak polarized emission is also detected 
in the outer jet, at the northern edge, with x ~ 0°. This 
change of the electric vector, roughly parallel to the jet 
direction at D and G7 contrasts with its value roughly 
perpendicular to it at the edge of the outer jet. The spec- 
tral analyses have shown that the jet at distances >1.5 mas 
of the core is an optically thin region, where also a north- 
south gradient in a was present. 

At 5 GHz, the distribution of the polarized emission 
differs significantly from that seen at higher frequencies. In 
Fig. |2|, we present the image obtained for the 3rd epoch 
of our monitoring. The two earlier epochs display very 
similar features in the jet. The core is less polarized (m 
1%) than at higher frequencies, with its electric vector 
obhque to the jet direction. Both facts may be caused by 
the blending of two or more emission regions with different 
polarization angles. The degree of polarization at 5 GHz in 
the jet reaches ~15%. The predominance of the polarized 
emission at the jet edges (especially at the northern one) 
is more obvious here than at 8.4 GHz. 



These findings seem to be consistent with the conclu- 



sions of Brown et al. (1994), who reported, at 5 GHz, a 
weakly polarized core in 3C 345 and a fractional polariza- 
tion reaching 15% in the jet. The electric vector position 
angles were found to be variable in the core from one epoch 
to the other but perpendicular to the jet for all epochs. 
VLBA observations of 3C 345 at 8.4 GHz made in 1997.07 



(Taylor 1998) yield an image that is very similar to our 
image at 1996.81, and both images imply a magnetic field 
aligned with the jet direction at 3-10 mas distances from 
the core. Cawthorne et al. ( 1993| ) suggested that the longi- 
tudinal component of the magnetic field can increase with 
the distance from the core as a result of shear from the 
dense emission line gas near the nucleus. At larger dis- 
tances from the core the shocks may become too weak to 
dominate the emission, resulting in the observed electric 
field perpendicular to the jet. This is possibly the physical 
mechanism that can explain the phenomena observed in 
our maps. 

We also investigated the possible existence of a differ- 
ential rotation measure (RM) at the edges of the outer 
jet, comparing cuts across the jet at 5 GHz and 8.4 GHz. 
Such a RM would be caused by a toroidal magnetic field 



(shear-like) in the jet (see Gabuzda 1999). Owing to the 
uncertainties in x, the only reliable result is that the values 
of the EVPA are ^ 40° bigger in the southern than in the 
northern edge, both at 5 and at 8.4 GHz. No conclusions 
about the presence of this RM can be extracted. 

4. Summary 

We have observed the supcrluminal QSO 3C 345 with 
VLBI at three epochs within one year, observing with 
the VLBA at four frequencies. Our results quantify the 
supcrluminal motions of components C8 and G7 with re- 
spect to the core component D, and they show a remark- 
ably complex polarization structure near the core, which 
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Fig. 12. Polarized intensity 
electric vectors (x, length pro- 
portional to p, 1 mas in the map 
is equivalent to lOmJy/beam) 
overlaid on total intensity (/) 
contours (3mJy/beamx — 
1,1,2.24,5,11.18,25,...) and 
grey scale polarized intensity 
(p, grey scale up to the peak 
of brightness, 40.5 mJy/beam) 
images for 3C 345 at 5 GHz, 
epoch 1996.81. It is obvious 
that the electric vector is almost 
perpendicular to the jet at core 
separations from 3 to 7 mas. 
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Fig. 13. Summary of physical 
properties observed in the two re- 
gions of the parsec-scale jet of 
3C345. The values for a and m 
refer to our four observing fre- 
quencies (5, 8.4, 15, and 22 GHz), 
in general. 



provides evidence for emerging components and changing 
projected jet direction within 3 mas from the core. A com- 
parison of our total intensity data at the four frequencies 
with previous model fitting results will provide a better 
understanding of the kinematical properties of the jet. 

Spatial blending may be partly responsible for the low 
core polarization observed for quas ars in centimeter VLBI 
polarimetry. Roberts et al. ( |l987 ) mention the possible 
interference between components with different x i^i the 
polarized intensity for the case of OJ 287 at 5 GHz. Similar 
features can be present for the case of 3C345, especially 
at 5 GHz, since different components have different values 
of X at higher frequencies. 



The polarized intensity images show alignment be- 
tween the jet direction and the electric vector position an- 
gle in the inner regions of the jet at high frequencies, and 
orthogonality of the electric vector and the jet direction in 
the outer regions of the jet. These results are fully consis- 
tent with the previously reported observations of Brown et 
al. ( |1994| ), Leppanen et al. ( |1995D , Taylor ( |199^ ), covering 
different frequencies and jet regions in 3C345. 

The interpretation of Lobanov & Zensus ( 1996| ) about 
a transition in the jet can be confirmed with the scenario 
provided by our observations (see Fig. |l3| ). The core re- 
gion (up to 1.5 mas), which presents magnetic field orien- 
tations predominantly transverse (Figs. ^ & |ll|), would be 
dominated by strong shocks (as components C8 and C7). 
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Fig. 11. VLBA /, p and x images of 3C345 at 15 
and 8.4 GHz, epoch 1995.84. The total intensity / is 
represented with contours (value of 6mJy/beamx — 
1,1,2.24,5,11.18,25,...), over a grey scale polarized 
intensity map (peaks of brightness of 169.3 and 
141.4mJy/beam, respectively, grey scale normalized to 
the common maximum polarized flux density -at 15 GHz- 
), and the electric vectors (x, length proportional to p, 
1 mas in the map is equivalent to 50mJy/beam). 

This shock dominance is also consistent with the peaked 
spectrum (Fig. 0). At larger distances (beyond 1.5 mas), 
the trajectories of the components are less curved than in 
the inner jet and show faster speeds, the magnetic field is 
longitudinal (Fig. |l^) and the spectrum is smoother (Fig. 

0)- 

The jet in 3C 345, thus, apparently manifests different 
behaviors of the magnetic field, which can be explained in 
terms of different mechanisms: (i) magnetic field compres- 
sion in the shocked, inner jet, (ii) ordering of the magnetic 
field in a shear layer and edge-brightened polarized emis- 
sion in the outer jet, (iii) decoupling of two fluids within 
the jet with different properties, obliqueness of shocks in- 
teracting with the outer edge of the jet (C8), etc. To disen- 
tangle the different physical phenomena that are present 
in the jet, further observations are needed, with the goal of 
higher resolution, and including polarimetry. The VSOP 
facility and the use of mm-VLBI will provide important 
clues to the innermost structures of the radio source (Klare 



et al. in preparation). Continuing the multi-band monitor- 
ing of 3C 345 with VLBI, enhanced by regular polarimet- 
ric studies, and also observations at higher frequencies and 
with better resolution (like those provided by orbital VLBI 
results with HALCA) should help to better constrain the 
models of this QSO. 
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